Critical exponents for random knots 
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The size of a zero thickness (no excluded volume) polymer ring is shown to scale with chain length N in the same 
way as the size of the excluded volume (self-avoiding) linear polymer, as AT", where v « 0.588. The consequences 
of that fact are examined, including sizes of trivial and non-trivial knots. 
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Knots are encountered in virtually every branch of 
physics, as evidenced by the series of books on "Knots 
and Everything." Knots keep entertaining physicists 
and matematicians for one and a half century, since 
W.Thomson [uj. In a more specific context, knotted 
polymers and DNA remain in the center of attention for 
over three decades [§-Q. Nevertheless, many simple ba- 
sic questions remain unanswered. For instance, the most 
fundametal physical property of any macromolecule - its 
average size R(N), scaling with chain length N, and de- 
pending on the solvent conditions - is not understood for 
knotted polymers. In particular, the value of gyration 
radius determines diffusion, light scatterring, gel pene- 
tration and several other properties of macromolcculcs. 
It is known in details for a variety of linear polymers un- 
der different solvent conditions. The mere closing ends 
and making a ring-shaped chain instead of a linear one re- 
duces its size only marginally, by a factor of y/2 for Gaus- 
sian chain and certainly preserving the power law even for 
excluded volume one. The really strong dramatic effect, 
as we show below in this paper, is due to the quenched 
restricted set of conformations - those which are topo- 
logically equivalent. Indeed, topological type must be 
quenched for all real closed ring-shaped polymers whose 
segments are unable to pass through each other. The 
only exception is DNA in the presence of topo-II en- 
zyme B, and the very fact that Nature has invented 
such marveluos tool indicates how serious is the topolog- 
ical problem for long polymers. Thus, the formulation 
of the problem is very simple: given the polymer ring 
of a quenched topology, we want to estimate its average 
size. So far, there have been only very few works on this 
problem P,R|, and overall the problem remains poorly 
understood. The purpose of this note is to propose some 
heuristic arguments addressing this problem. 

1. In this paper, we consider the simplest ring poly- 
mer with no excluded volume ||, consisting of some N 
freely-jointed straight segments, length £ each. Assume 
that every spatial shape, or conformation, of this polymer 
is just as likely as any other. In mathematical language, 
our "polymer" is simply a closed broken line embedded 
in 3D. With probability 1, this object does not have self- 
intersections, and, therefore, represents a knot. It may be 
a trivial knot, if it can be continuously and without self- 



crossings transformed into a plane polygon topologically 
equivalent to a circle, or it may be one of the non-trivial 
knots. It has beem proven as early as in 1988 p(J that 
the probability of a trivial knot configuration decays ex- 
ponentially with the chain length, N: 



Ptk(iV) ~ exp(-JV/JV ) 



(1) 



where "tk" subscript stays for "trivial knot" and a con- 
stant Nq represents the characteristic length. This math- 
ematical prediction has been beautifully confirmed in 
computer simulation (jf!), yielding No ~ 340 ± 4. 

The restriction to broken-line polymer of straight seg- 
ments is not necessary. Alternatively, one can think, 
for instance, of a smooth worm-like ring object of the 
length N£ with effective segment I [n2[. One can also 
make straight segments with Gaussian distributed ran- 
dom lengths pa]. The exponential result (nh holds for 
these models as well as for all similar ones, except per- 
haps with different values of No- However, one cannot 
use continuous model (the Wiener trajectory), since it 
corresponds to the limit N — * oo, £ — * 0, N£ — const, 
where essentially all conformations are heavily knotted, 
or trivial knot probability is exactly zero. 

While rigorous proof of the result (|l]) is not simple, it 
can be readily understood on a hand-waiving level. In- 
deed, consider our polymer as a sequence of N/g blobs, 
with g segments in each blob. The probability of trivial 
topology for each blob is Ptk(sO> and formula ([[]) is equiv- 
alent to saying that ptk(-A) = [ptk(ff)] > which means 
that topological constraints between blobs are negligi- 
ble compared to that between segments inside each blob. 
The temptation is then to conclude that vast majority of 
knots occur on the small length scales. However appeal- 
ing, this jargon statement is confusing: indeed, while blob 
may be much smaller than the entire polymer, g <C N, 
the argument above only holds if it is large compared to 
No (<? ^> No), which itself is numerically large. 

Exponentially small probability of a trivial knot (Fy) 
immeaditely explains why topological constraints can sig- 
nificantly alter the polymer size. Indeed, £N X I 2 is the size 
averaged over all conformations. However, the size of a 
trivial knot is an average over an exponentially small sub- 
set of conformations. When relatively compact complex 
knots are removed, the remaining average increases. In 
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order to bring this idea to a quantitative level, a couple 
of further preliminary arguments will be handy. 

2. Consider now two ring polymers, with iV segments 
each. They form a link, trivial or non-trivial If rings 
overlap, i.e. the distance between their centers of mass 
is smaller than their sizes, the probability that the link 
is trivial vanishes, pt\{N) — > 0, when N — » oo. While the 
exact asymptotic expressions for pti(N) are not known, 
it should be qualitatively similar to the trivial knot prob- 
ability (Fy), and we assume that the corresponding char- 
acteristic chain length iVi is of the order of N . 

Speaking about the trivial link probability, one has to 
make the distinction between a link made by two rings 
each of which may be an arbitrary knot, and a link made 
by two trivially knotted rings. Since the former probabil- 
ity is higher, and goes to for large N, the latter, which 
we use below, goes to zero as well. 

3. Since trivial link topology is highly unlikely for 
overlapping rings, the untangled ring polymers effectively 
present excluded volume for each other, even if there is 
no bare excluded volume for monomers. Simple scaling 
argument suggests that the excluded volume can be es- 
timated as 



v cxcl (N)~R 3 (N){l~ Ptl (N)} 



(2) 



where R{N) is the (properly averaged) size of one ring. 
The concept of "topological excluded volume" was intro- 
duced in one of the first works on computer simulation of 
knots and links |L5| . It can be easily understood in terms 
of an exactly solved model fl(||l7j in which one polymer is 
a ring and the other is a straight line. It has recently been 
rigorously proven |L8[ that the excluded volume scales as 
R 3 (N) in N — > oo limit for two rings if their Gaussian 
linking number remains zero, even though this condition 
does not guarantee the untangled topology. As regards 
eq. (0), it is actually almost trivially correct. Indeed, 
1 —pti(N) is the probability for two phantom rings to be 
tangled, i.e., to adopt a conformation prohibited for real 
(non-phantom) untangled rings 119]. 

4. We are now prepared to directly attack the question 
of the size, Rtk(N), of a trivial knot without excluded vol- 
ume. Let us color two different pieces of our iV-segment 
ring, both containing some g <C N monomers (segments). 
Since neither of the pieces has open ends, they present 
uncrossable objects for one another. Therefore, the ar- 
guments above apply and suggest that these two pieces 
exclude for each other some volume which is of the or- 
der of v exc \(g). While this identification of the pieces of 
a ring chain with separate small rings may be doubtful 
for small g, we expect it to become increasingly accurate 
with increasing g, especially when g becomes greater than 
No- In that case, pti(g) is negligible, and we end up with 
^exci(ff) — R 3 (g) for g > Nq, which immediately implies 
that the chain of g-blobs belongs to the universality class 
of self- avoiding walks, and thus 

JZt k (iV) ~ R(g) {N/gf , where v w 0.588 « 3/5 . (3) 



We expect this to hold at g > Nq, while at smaller scales 
the effect of topological constraints is only marginal. 
Thus, we assume R(g) ~ ^S 1 ^ 2 f° r 9 U P to about iV 
and apply the above result (|3|) choosing g ~ Nq: 



Rt*{N) 



£N r 



IN 1 / 2 if N <N 



N v if N > N 



-"+1/2 A ™ 



(4) 



Thus, chain size in the case of trivial knot topology is con- 
trolled by the standard excluded volume exponent close 
to 3/5. This confirms the conjecture made by J. des 
Cloizeaux as early as in 1981 M. We see also that this 
"swollen" regime comes only for really long chains, be- 
cause No is numerically large jflj. Furthermore, while 
in terms of parameters the cross-over to the "swollen" 
asymptotics is controlled by No, there are grounds to ex- 
pect that this asymptotics actually develops at N around 
a few times No (which we expressed with the 3> sign in 
formula (Q)). On the other hand, although topological 
constraints lead to swelling only when N is really large, 
the prefactor in the second line of formula fl) is of order 
unity: A^ +1/2 » 34(r a088 « 0.6. 

It is an interesting challenge to produce a renormal- 
ization group analyzis for the trivial knot polymer. One 
of the difficulties is that regular e-expansion techniques 
are hardly applicable, since topological constraints do not 
exist for linear polymers in dimcsnions other than 3; as 
a matter of fact, these constraints exist in the space of 
dimension d for the objects of dimension d— 2. The man- 
ifestation of that difficulty is seen in the fact that trivial 
link probability enters the expression (0) for the scale- 
dependent excluded volume. On the other hand, such 
RG calculation, once achieved, can shed a long awaited 
light on why Nq is as large numerically as it is. 

5. Consider now some non-trivial knot, still without 
any excluded volume. This can be addressed using the 
so-called "tube inflation" or "ideal knot" approach. This 
idea has been proposed in several contexts by different 
authors, in the aspect addressing knot entropy it was 
suggested in [£0|. Recently, the approach received a sig- 
nificant attention |2l| . In particular, the present author's 
contribution in this book addresed the size of an arbi- 
trary knot. However, the result (||) was not established 
at the time; instead, the size of an arbitrary knot was es- 
timated under the assumption that size of a trivial knot 
scales as N** with an unknown exponent /x. We can now 
directly use the results of the author's work in pi] , sim- 
ply substituting fi — v there. To make the present work 
self-contained, we repeat the argument briefly. 

When polymer is a knotted ring, fluctuations of each 
segment are restricted by the neighboring segments. Fol- 
lowing the standard trick of the reptation theory J2^] , we 
replace these constraints by an effective tube confining 
the entire polymer. The difficult part of the problem is 
that, unlike the usual case of a melt p2j |, the tube shape 
for the swollen single chain fluctuates very strongly, and 
we have to find a reasonable approximation for it. To this 
end, we argue that there are two possibilities to achieve 
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maximal entropy, corresponding to "uniform" and "phase 
segregated" structures. Let us examine both. 

We call the structure "uniform" when the freedom 
of transverse fluctuations is about the same everywhere 
along the polymer. That corresponds to finding the max- 
imal diameter, D, for which the tube of the length L can 
be knotted into the given knot without self-penetration. 
The shape of this "maximally inflated" tube, or the shape 
of its central axis line, is called "ideal representation of a 
knot." It is characteristic of the knots topology, and so 
is this tube "axis ratio" p = L/D. 

Thus, we replace the real topological constraints with 
the confinement within the tube whose aspect ratio, p, 
is that characteristic for the given knot topology. Tube 
is closed, and polymer makes exactly one turn within 
the tube. We assume further that the polymer is not 
knotted inside the tube p3fl . These two conditions to- 
gether guarantee that the polymer has exactly the desir- 
able topology. We can now determine overall size of the 
polymer in space, R, using the following Flory-style argu- 
ment: When R increases, polymer loses entropy because 
it gets stretched along the tube axis; when R decreases, 
polymer loses entropy because it gets squeezed perpen- 
dicular to the tube axis. Thus, the equilibrium size R 
can be conjectured to balance these two entropic factors. 
To implement this idea, we choose D and L such that 
LD 2 ~ R 3 and, since L/D = p, that means L ~ Rp 2 / 3 , 
D ~ Rp~~ x / 3 . We further recall that the unknotted poly- 
mer behaves just as an excluded volume one on large 
scales. That allows us to employ the standard scaling 
arguments Q to estimate the relevant entropy. Specifi- 
cally, transverse compression and longitudinal stretching 
in the tube are associated with, respectively, concentra- 
tion blobs and Pincus tension blobs: 



S~- 



L 

iV 7 - 



V(i-") 



/ N 3v 

Vld 2 



1/(3^-1) 



(5) 



Provided the above expressions for L and D, this entropy 
has a maximum with respect to R at R ~ EN"p~ v+1 / 3 , 
and we identify it as the average size of the knot. This re- 
sult applies as long as each of the aforementioned blobs 
remains larger then Nq monomers. Equivalently, equi- 
librium tube diameter D should be greater than £Nq , 
yielding N > pNq. For smaller blobs, Gaussian statistics 
should be valid, which means physically that we can con- 
sider a phantom |19] polymer confined within the tube. 
This case has been already considered in [g0| . Collecting 
everything together, we obtain 



R k (N) 



Wr 



-"+1/2 



Ofi/ip-V* if N < Pk N 
N»pl 



Pk 

"+1/3 



if N > Pk N 



> (6) 



where we have included the subscript in p k to emphasize 
that this parameter is taken for the given knot k. 

Thus, in terms of A-dependence in the limit of very 
large N, every knot acts as an effective excluded volume, 
making the exponent equal to v w 0.588. However, more 



complex knots, with larger p, are significantly less ex- 
panded than trivial or simple knots. Furthermore, for 
complex knots the onset of swollen behavior is pushed to 
the range of really long chains. These results are consis- 
tent with numerical observation |0,Gq] suggesting that 
probability of any particular knot, and not only trivial 
one, decays exponentially at large N (|l|), with the char- 
acteristic length Nq independent of the knot type. Qual- 
itatively, one can say that increasing N for the polymer 
with any given topology eventually leads to the situa- 
tion when polymer is dominated by a very long strings in 
which knot restriction is not felt at all, and thus it should 
become in this sense equivalent to the trivial knot. 

Another way to look at the eq. (g) is from the point of 
view of p dependence. For the very long chain, with N 3> 
./Vo, topological "memory" acts as an excluded volume 
for relatively simple knots, with p < N/Nq, for which 
R obeys the last line of eq. (W). With increasing knot 
complexity, the system crosses over to the other regime, 
described by the upper line in eq. (n). Note that at the 
cross-over the chain size is about R ~ e.N 1 l 3 N G ' . In 
terms of TV-dependence, this is already a collapsed chain. 

Note also that the maximal value of entropy (0) is 
about S ~ — p: the knot creates about one constraint 
per polymer length which is about tube diameter. 

6. As we mentioned, there is an alternative way to 
maximize entropy, which we call "phase segregated." To 
explain what we mean by that, consider first a prime 
knot. It may be entropically favorable to tighten this 
knot as much as possible, losing virtually all the free- 
dom for segments involved in the tight knot, but gaining 
maximal entropy for the rest of the chain which then fluc- 
tuates as a free unknotted loop. For the composite knot, 
this picture must be slightly modified: we still imagine a 
long free unknotted loop with several little knots - prime 
componets of the original composite knot - independently 
tightened in different places on the loop. 

To estimate roughly the entropy of such knot segre- 
gated state, we can use again formula (||), applying it 
twice (or several times for a composite knot): for the 
tightened part, where p is roughly equal to that of the 
original knot, and for the unknotted loop part, where p is 
of order one. The chain lengths for these two are about 
p and N — p, respectively. Since equilibrium entropy is 
linear in p, this simplest estimate yields that phase segre- 
gated state is about as favorable entropically as the uni- 
form one. A more sophisticated approach may be needed 
to decide which of them is more stable. It is also likely 
that the answer is sensitive to the details of local chain 
geometry (eg, freely-jointed vs worm-like segments, and 
the like). 

7. Although the qualitative statement that pti{N) is 
practically zero for long rings is sufficient for the main 
stream of arguments in this paper, to estimate pti(N) is 
itself an interesting problem. The simple estimate is ob- 
tained in the following way. Consider two strongly over- 
lapping rings. Since we are counting probabilities over all 
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possible configurations of rings, regardless of their knot- 
tedness, the dominant statistics is Gaussian, and there 
should be about ^/N contacts between the rings. As 
long as A is not very large, each contact can be thought 
of as capable to provide up to one (positive or negative) 
unit of linking. Therefore, there are about VA differ- 
ent linking states, with the probability of order 1/yN 
for each of them, including the untangled one. When A 
gets really very large, the situation becomes more compli- 
cated, because each of the vN "contacts" is itself a large 
crowd of segments, which can realize a variety of linking 
arrangements. Since entanglements in different contact 
regions do not commute to each other in this regime p4[ , 
there appears exponential diversity of different types of 
linkings. Adjusting prefactors for continuity, we arrive at 



Pti(JV) 



exp 1 - ^/N/N- 



^/NtJN if N < Ai 
if N > Ax 



(7) 



where subscript "tl" stays for the trivial link, Ai is the 
characteristic length. The first line of this formula is in a 
very good agreement with the simulation data of the work 
p5| where links of the lengths up to 500 were examined. 
This is consistent with our conjecture that Ai is of the 
order of Ao. It would be interesting to simulate longer 
chains, as we expect the iV-dependence to cross-over to 
exponential in the range of several hundered segments. 

It is also interesting to present the blob-style argument 
to explain the lower line of the eq. (M) . Consider each ring 
as consisting of N/g blobs, with g segments in each blob. 
About \f N / g of the blobs from one ring are expected to 
overlap with the blobs of the other ring. Now the bold 
assumption comes (supported by an admittedly vague 
argument): in order for two rings not to be linked, all 
of the blobs must be untangled. Since the "dangerous" 
blobs are only the overalpping ones, we obtain p t \(N) = 

[Pti(d)} i which yields the familiar answer (Q). 

When this work has already been written, the paper pi] 
appeared. It reports simulation data for the ring-shaped 
polymer which preserves the topology of a trivial knot 
and does not have any excluded volume. The average 
size of this polymer has been measured for a series of 
lengths A = 2™, ranging from 32 to 2048 monomers. 
The data are consistent with our results. Specifically, the 
sizes of A-mers were found to follow Gaussian statistics, 
A 05 , for small A (n = 5,6,7), then power was slowly 
increasing with A, and finally the ratio of the sizes of 
2048- and 1024- mers was close to 2 - 58 . This seems in 
perfect agreement with the result eq. (|J). 

The author has a pleasure to acknowledge useful dis- 
cussions with B.Duplantier, I.Erukhimovich, S.Nechaev, 
D.Nelson, TWitten. 
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